Realization of a spin-liquid ground state in the high-pressure phase of Ba3YIr209 



Tusharkanti Dey,i S. K. Panda,^ A.V. Mahajan,i'EI I- Dasgupta,^'^ 
B. Koteswararao."* F. C. Chou,^ A. A. Omrani,^ and H. M. R0nnow^ 

^Department of Physics, Indian Institute of Technology Bombay, Powai, Mumbai 400076, India 
^ Centre for Advanced Materials, Indian Association for the Cultivation of Science, Jadavpur, Kolkata-700032, India 

^Department of Solid State Physics, Indian Association for 
the Cultivation of Science, Jadavpur, Kolkata 700032, India 
'^Center for Condensed Matter Sciences, National Taiwan University, Taipei 10617, Taiwan 
^ Laboratory for Quantum Magnetism, Ecole Polytechnique Federale de Lausanne (EPFL), CH 1015, Switzerland 

Abstract 

Ba3YIr209 crystallizes in a hexagonal structure and transforms into a cubic structure under an applied pressure of 8GPa 
at 1273 K. While the ambient pressure (AP) sample undergoes long-range magnetic ordering at ~ 4K, the high-pressure (HP) 
synthesized sample does not order down to 2K as evidenced from our susceptibility, heat capacity and nuclear magnetic 
resonance (NMR) measurements. Further, for the HP sample, heat capacity measurements demonstrate the presence of a 
linear term in C/T vs. in the range of 2 - 10 K with 7 = lOmJ/mol-Ir K^. The NMR shift has no temperature (T) 
dependence in the range of 4 — 120 K and spin-lattice relaxation rate varies linearly with T in the range of 8 — 45 K (above which 
it is T- independent). Resistance measurements confirm that both the samples are semiconducting. Our ab-initio calculations 
reveal that strong spin orbit coupling (SOC) plays a crucial role in determining the magnetic and insulating properties of this 
system in both the phases. We suggest that the cubic (HP) phase of BasYIraOg is a 5d based, gapless, spin-liquid. 



PACS numbers: 76.60.-k, 71.15.Mb, 75.70.Tj 

The 3d transition metal oxides have been studied ex- 
tensively to explore novel properties like superconduc- 

[ tivity [3, 01 colossal magnetoresistance etc. In these 
materials, the orbital angular momentum of the valence 

] electrons is largely quenched and a large onsite Coulomb 
energy often drives the materials to a Mott insulating 
state. In contrast, for 5d group based materials, this on- 
site Coulomb energy is smaller by an order of magnitude 
and one could expect more metallic and less magnetic 
behavior. However, some 5d Ir-based materials such as 

■ Na4lr308 0, Sr2lr04 Si], Na2lr03 etc. are found to 

[ be insulators showing exotic magnetic properties. These 
materials are insulating due to the presence of strong 
spin-orbit interactions and are often called as spin-orbit 
driven Mott insulators. Such materials are relatively less 
explored and expected to offer a promising playground 
for materials researchers. 

To explore this spin-orbit driven physics further, we 
have recently investigated hexagonal Ba3lrTi209 and 
suggested it to be a spin-orbit driven liquid based on 

[ a triangular lattice Ba3YIr20g has a similar chemi- 
cal formula and it crystallizes in the hexagonal structure 
(P63/mmc) with Ir-Ir structural dimers arranged in an 
edge shared triangular fashion Q. Since all the Ir are 
equivalent, they should have a fractional oxidation state 
of -|-4.5 in a simple ionic picture. This should presum- 
ably lead to metallic behavior in sharp contrast to the 
insulating behavior seen in resistivity measurements. In 
spite of the triangular lattice being frustrated, Doi et al. 
[Sl reported magnetic ordering in the above compound at 
~ 4K. It is likely that there is significant spin-orbit cou- 
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pling in Ba3YIr209 which therefore warrants a deeper 
investigation for possible evidence of an exotic ground 
state. In a similar system Ba3NiSb209, the ambient pres- 
sure sample orders at 13.5 K but its structure changes 
when exposed under high pressure and it has been pro- 
posed to become a spin liquid It would therefore 
be interesting to explore Ba3YIr2 09 under high pressure 
to see if any change in crystal structure and hence in 
magnetic properties takes place. 

Herein we report x-ray diffracion, magnetization, heat 
capacity and ^^Y NMR measurements accompanied 
by first-principles electronic structure calculations for 
both ambient (AP) and high (HP) pressure phases of 
Ba3YIr209 with focus on the HP phase. The AP sample 
is hexagonal, insulating and orders at 4 K as reported 
earlier. On the other hand, the HP sample has a cubic 
structure, exhibits a Curie susceptibility but does not 
order down to 2K. More specifically, the HP sample is 
semiconducting but has a metal-like linear heat capacity 
coefficient 7 = lOmJ/mol-Ir K^. The ^^Y NMR shift is 
independent of temperature T below 120 K and the ^^Y 
NMR spin-lattice relaxation rate crosses over from T- 
independent behavior at high temperature to a linear T- 
dependence below about 45 K. These results point to the 
existence of low-energy excitations at low-temperatures. 
In the absence of metallic behavior in the resistivity, the 
low-T data suggest the formation of an exotic ground 
state, possibly a gapless spin-liquid state. The electronic 
structure calculations reveal the important role of spin- 
orbit coupling for the description of the insulating state 
in both the phases. Together these results establish a 
new spin-orbit driven Mott-Hubbard insulator. 

We have prepared polycrystalline sample of Ba3YIr209 
by solid state reaction method as detailed in Ref. 
Further, the AP sample was treated under 8 GPa pressure 
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Figure 1: X-ray diffraction pattern of the HP phase of 
Ba3YIr209 sample is shown along with a two phase refine- 
ment with space group Fm-3m. The (hkl) positions are also 
indicated. A 5% impurity phase of ureacted Ir is found. The 
peaks corresponding to Ir are marked with (*). The main 
peak of the primary phase is shown in the inset. 

at 1273 K for 30min to get the HP sample. Due to the 
small size of the HP cell, only about 150 mg of sample 
was obtained. 

Our XRD analysis confirmed that the AP sample was 
formed in a hexagonal structure (PGs/mmc) as reported 
earlier Q but after high-pressure treatment the struc- 
ture changed to a cubic double perovskite (Ba2MM'06) 
phase with space group Fm-3m. It is relevant to men- 
tion that Ba3NiSb2 09 was also reported to transform 
from a hexagonal (PGa/mmc) to cubic (Fm-3m) phase 
when treated under 9 GPa pressure at 1273 K for 30 min 
[13. In the XRD pattern (Fig. [J) of the HP phase, 
a few extra peaks (marked as *) are present which are 
identified as due to unreacted Ir and account for about 
5% of the contribution. Also, all the major peaks have 
a shoulder on the right side. These shoulders possibly 
arise from a double perovskite phase with a slightly dif- 
ferent cell parameter. Accordingly, a two-phase Rietveld 
refinement (Fig. [T]) performed with our XRD data yields 
lattice parameter a = 8.3248A for the primary phase 
and a' = 8.202A for the secondary phase (about 10% 
content). Refinement with a cubic perovskite structure 
(chemical formula Ba(Y]^/3lr2/3)03) with approximately 
half the lattice parameter of the double perovskite (chem- 
ical formula Ba2lr(Y2/3lri/3)06) did not give satisfactory 
results. We find that an ordered arrangement of IrOg and 
Y2/3lrj/306 octahedra form the cubic double perovskite 
Ba2l r(Y2/3 lri/3)Og structure of the HP phase (shown in 
Fig. [Mil] of Supplemental Material (SM) We note 

that every Y will have 6 Ir as its nearest neighbors while 
as next nearest neighbors it will have 33% Ir and 67% Y. 

Inset of Fig. [2]shows the magnetic susceptibility of the 



AP as well as HP samples measured as a function of T in 
a field 5 kOe. For the AP sample, magnetic susceptibility 
shows a kink at ^ 4K indicative of a long range mag- 
netic ordering as reported earlier Q. Our calculations 
based on GGA+U+SOC, pertaining to the AP sample, 
suggest that among the calculated structures, the one 
with ferromagnetic intra-dimer coupling and an antifer- 
romagnetic inter-dimer coupling (AFM2) has the lowest 
energy (see SM [llj ) . Susceptibility data of the AP sam- 
ple could be fitted to the Curie law in the temperature 
range 30 — 300 K which yields a T-independent suscepti- 
bility xo = 3.8 X 10~^cm^/mol Ir and a Curie constant 
C = 0.0125 cm^K/mol Ir. Below about 30 K, the suscep- 
tibility deviates from Curie behavior due to building up 
of magnetic correlations. For the HP sample, no anomaly 
indicative of a magnetic ordering was seen down to 2 K 
nor was there any ZFC-FC splitting observed. Data for 
the HP sample could be fitted to the CW law in almost 
the whole temperature range resulting in the parameters 
Xo = 3.9 X 10-'*cm3/mol Ir, C = 0.0045 cm^K/mol Ir 
and 9 ~ — 1.6K. The small 9 found for the HP sample 
in the full temperature range indicates the weakness of 
magnetic interactions. The C value for the AP sample 
is nearly a factor of 30 smaller than for isolated 5* = 1/2 
magnetic moments (0.375 cm'^/mol). For the HP sample, 
the C value is smaller by a further factor of 3 compared 
to its AP counterpart while the xo is nearly unchanged. 
The absence of long-range order in the HP structure is 
also anticipated by our electronic structure calculations 
which find various magnetic structures to be degenerate 
(see SM [ll|)- Considering xo = Xcore + Xw, where Xcore 
for our sample composition is —1.37 x 10~^cm^/mol Ir, 
we get Xvv = 5.3 x 10"'' cm^/mol Ir. This is nearly within 
the range (1 x 10~^ — 5x 10~^ cm^/mol Ir) of Xw observed 
in different Ir-based insulating oxides 0,111. Absence of 
a XPauii term suggests the sample is insulating. 

The low-temperature heat capacity can often throw 
light on the presence of low-energy excitations in the 
material. If one is able to separate out the phonon con- 
tribution, the residual contribution should represent the 
electronic part which is generally of interest. Whereas in 
the AP sample a clear anomaly at about 4K is seen in 
the heat capacity (see Fig. [5]) which reflects long-range 
order, the anomaly goes away in the HP sample. This 
indicates that magnetic ordering has been suppressed by 
the application of pressure. It is further interesting to 
note that for the HP sample below about 10 K, Cp/T vs 
is linear (see Fig. ^ and gives a non-zero Sommerfeld 
coefficient 7=10 mJ / mole-Ir and (3 = 0.44 mJ / mol-Ir 
K**. Although this 7 value is typical of metallic materi- 
als, our resistance measurement on single crystals show 
semiconductor like behavior (see SM [111). These results 
could imply the exotic possibility of the absence of charge 
excitations while spin-excitations are present. In such a 
scenario, one can interpret the linear contribution to the 
low-T heat capacity as coming from excitations from a 
spinon Fermi surface as might happen in a gapless spin- 
liquid. 
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Figure 2: Gp/T vs is shown for the AP and HP samples. 
The solid line is a linear fit to data for the HP sample. Inset: 
Magnetic susceptibility of AP and HP samples measured in a 
5kOe field and their fitting (solid lines) with Curie-Weiss law. 
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Figure 3: PDOS of Ir-d states within GGA+U+SOC for (a) 
AP and (b) HP phase. Schematic to explain the insulating 
behavior of AP phase is shown in the inset. 



Before we present further evidence of such non- 
standard behavior of the HP phase in other measure- 
ments we shall discuss the origin of the unusual insulat- 
ing state in both the AP as well as HP samples where Ir 
is expected to have a fractional oxidation state of -1-4.5. 
To understand this issue we have performed electronic 
structure calculations within local (spin) density approx- 
imation with generalized gradient correction (GGA) of 
Perdew-Burke-Ernzerhof including Hubbard U and spin 
orbit coupling (SOC) (for details see SM [uj). In Fig 
[3Ja) we have displayed the partial density of states of 
Ir-d states in the AP phase for the AFM2 structure. 
The DOS reveal that the system is insulating. In the 
following we shall argue that in contrast to other irri- 
dates, in addition to SOC and Hubbard U the strong 
intra-dimer coupling plays a key role for the system to 
be insulating. The DOS reveal that the exchange split- 
ting is much smaller compared to the crystal field split- 
ting. So Ir [d^-^) remains in the low spin state. As a 
consequence majority t^g states are completely occupied 
and the remaining 1.5 electrons/Ir are accommodated to 
the minority spin channel. The origin of the gap in the 
minority spin channel is schematically illustrated in the 
inset of Fig ISJa) where we have only shown the Ir-t2g 
states in the minority spin channel for an Ir-Ir dimer. 
The threefold degeneracy of the i2g states is lifted by 
the distortion of IrOe octahedron and by SOC (see SM 
[m for details). This is however not enough to open up 
a gap in the minority spin channel with 1.5 electron/Ir. 
Interestingly a strong intra-dimer coupling produces a set 
of bonding and anti-bonding states shown in the middle 
panel of the inset of Fig. [DJa). The three electrons avail- 
able per dimer will occupy the states and can open up a 
gap at the Fermi level provided the bandwidth is smaller 
than the splitting due to dimerization. In such a situa- 
tion, a small but finite Hubbard U (U = 3 eV) promotes 
an insulating state by increasing the separation between 
the occupied and unoccupied states. In Fig. Ellb) we 
have shown the total non-spin polarized DOS for a model 
system Ba2lr(Y3/4lri/4)06 which is very close to the ex- 
perimentally synthesized Ba2lr(Y2/3lr4/3)06 obtained in 
the HP cubic phase. Our calculations reveal strong spin- 
orbit coupling and difference in the filling of the i2g levels 
for Ir and Iry (i.e. Ir at the Y site) are responsible for 
promoting the insulating state. 

Having established the insulating state for the AP and 
HP phases, we looked for signs of novel magnetic be- 
havior using the *^Y NMR local probe via spectra and 
spin-lattice relaxation rate measurements (which probe 
the low-energy excitations) on both the AP sample and 
the HP one (a mere 150 mg). Here, we focus on the data 
below about 120 K which is relevant to the theme of this 
paper. The spectrum of the AP sample shifts negatively 
with decreasing temperature (inset (a) of Fig. 2]), in line 
with the Curie-like susceptibility and a negative hyper- 
fine coupling constant. The linewidth too, increases with 
decreasing temperature and tends to diverge at ~ 4 K 
(inset (b) of Fig. 0]) as a signature of the onset of long- 
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Figure 4: Y NMR spectra at different temperatures for the 
HP sample are shown. The arrow indicates measurement tem- 
perature in increasing order. Inset: (a) and (b) show the Shift 
and FWHM, respectively, of the AP sample as well as the HP 
sample as functions of temperature. The arrow in the in- 
set (b) indicates that FWHM for the AP sample goes much 
higher at low temperature (~ 40kHz at 4K). 
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Figure 5: l/Tj, for the HP sample is shown on the left axis. 
1/Tl increases linearly with T in the range 8 — 45K. Above 
45 K, 1/Tl is almost independent of T. Solid lines are guides 
to the eye. On the right axis, l/Ti for the AP sample is 
shown. 



range order. A corresponding signature is also seen in 
the divergence of the ^^Y NMR 1/Ti for the AP sam- 
ple. In contrast, the spectra for the HP sample (shown 
in Fig. U]) are narrower and the width increases only 
slightly with decreasing temperature without any diver- 
gence. This would seem to be in tune with the much 
weaker Curie-term and the lack of magnetic order seen 
in the susceptibility. Somewhat surprisingly, though, the 



NMR shift for the HP sample is T-independent below 
120 K. The ^^Y NMR relaxation rate (Fig. O for the HP 
sample is T-independent at high-T as might be expected 
in a paramagnet insulator but crosses over to a linear 
variation below about 45 K (the rate corresponding to 
the longer Ti component {1/Tl) is plotted here; (see SM 
[nj for details). The constancy of the NMR shift for 
the HP sample coupled with the Curie-like susceptibility 
variation in the same temperature range suggest one of 
the two things, (i) a very weak hyperfine coupling results 
in an insensitivity to the small susceptibility variation 
(this is in accord with the long Ti) or (ii) the low-T in- 
crease of the susceptibility is of extrinsic origin and that 
the intrinsic spin susceptibility saturates as one decreases 
temperature and remains constant below 120K. 

Our susceptibility, resistivity data and electronic struc- 
ture calculations support the paramagnetic insulator pic- 
ture but once again, a non-zero 7 and crossover in l/Ti 
data run counter to this idea. On the other hand, the 
linear behavior in Cp/T vs with a non-zero inter- 
cept and the linear variation of l/Ti with T, inspite of 
insulating behavior point to a gapless spin liquid state. 
In this context it is relevant to point out NMR results 
on other spin liquid systems. For the 5d-based spin liq- 
uid system Na4lr308, ^^Na NMR by Takagi et al. 
find a negligible temperature dependence of shift below 
about 100 K since the susceptibility was saturated. In 
this case the ^'^Na 1 /Ti was constant above 200 K below 
which it followed a power law with power 1 in the range 
2 — 10 K and with power 2.5 in the temperature range 
100 — 200 K. They attributed this low-temperature linear 
behavior to low lying spin excitations. Itou et al. [T^ 
also found almost no temperature dependence of shift in 
their NMR data down to 18.7 mK on the organic 
system EtMe3Sb[Pd(dmit)2]2 with a spin liquid ground 
state while 1 /Ti was found to be constant in some range 
above 1 K and varied as below 1 K. In the organic spin 
liguid system k-(BEDT-TTF)2Cu2(CN)3, Shimizu et al. 
[la | found no change in ^H NMR spectral position from 
36.1 K down to 32 mK while 1/Ti followed a power law 
at low-temperature, from which the authors concluded 
that the ground state has gapless low lying spin excita- 
tions. Our NMR data are qualitatively similar to the 
cases listed above, providing further justification for the 
spin liquid state at low T . It is worth mentioning that 
another B site ordered double perovskite Ba2YMo06 was 
reported to have an exotic valence bond glas s state [l^ 
or a collective (gapped) spin-singlet state |l7| . 

In conclusion, a strong SOC plays a crucial role in the 
insulating behavior of both AP and HP samples. The 
magnetic behavior of the AP and HP phases are markedly 
different. We have also shown that the AP sample orders 
magnetically at ~ 4K. The ordering vanishes and the 
magnetism becomes very weak in the case of the HP sam- 
ple. For the HP sample, the linear term in heat capac- 
ity and linear temperature dependence of relaxation rate 
could be arising from the excitations out of the spinon 
Fermi surface of a spin liquid with very weak moments. 
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From our low-T data, we suggest that the HP phase of 
Ba3YIr209 is a Sd-based, gapless, spin hquid. 
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Supplemental material for "Realization of 
a spin-liquid ground state in the 
high-pressure phase of Ba3YIr209" 



I. DETAILS OF EXPERIMENTS 

Powder x-ray diffraction (XRD) measurements were 
performed at room temperature with Cu radiation 
(A = 1.54182A) in a PANalytical X'Pert PRO diffrac- 
tometer using Si for calibration. Analysis of XRD data 
for the HP sample reveals that the sample is transformed 
to a cubic double perovskite structure as shown in Fig. 
ISM- II The red octahedra represent IrOe and among the 
violet octahedra, 1/3 are IrOe and 2/3 are YOe- 
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Figure SM-1: The crystal structure of the high pressure syn- 
thesized cubic phase is shown. The blue and yellow atoms 
represent Ba and O, respectively. The red octahedra are IrOe 
octahedra while the violet octahedra contain Ir or Y at their 
centres. 

Magnetization measurements (shown in the main pa- 
per) were carried out in the temperature range 2 — 400 K 
and field range — 7 T using a Quantum Design SQUID 
VSM. Heat capacity measurements were performed us- 
ing the heat capacity attachment of a Quantum Design 
Physical Properties Measurement System (PPMS) in the 
temperature range 1.8 — 300 K and field range — 9T. 

We have measured resistance on small single crys- 
tals [600yLtmx 150/^mx lOOyLtm (approx.)] (from the same 
batch on which other measurements are done) for both 
AP and HP samples using the four probe technique. The 
temperature dependence of resistance for both the sam- 
ples is semiconductor- like (Fig. ISM-2[) . 

For ^^Y NMR measurements, a Tecmag pulse spec- 
trometer and a fixed magnetic field 93.954 kOe obtained 
inside a room-temperature-bore Varian superconducting 
magnet were used. Variable temperature was obtained 
with the help of an Oxford cryostat and accessories. The 
^^Y nucleus has spin 1 = 1/2 (100% natural abundance) 
and gyromagnetic ratio 7/27r = 2.08583 MHz/T. 
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Figure SM-2: Normalised resistance of both AP and HP sam- 
ple as a function of temperature. Ro is the resistance at room 
temperature and the values are mentioned in the figure. 
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Figure SM-3: Recovery of the Y longitudinal nuclear mag- 
netization after a saturating pulse is shown for the HP sample 
at a few temperatures. 



We have measured spin-lattice relaxation rate (l/Ii) 
using a standard saturation recovery method following 
a. '^/2 — t — {'^/2 — tt) pulse sequence in the temperature 
range 4 - 300 K for the AP sample and 8 - 120 K for the 
HP sample. For the AP sample, nuclear magnetisation 
recoveries are single exponential as expected for 1 = 1/2 
nucleus. Relaxation time (Ti) for the AP sample varies 
from about 10 ms to about 400 ms as a function of tem- 
perature. For the HP sample the nuclear magnetisation 
recovery after a saturating pulse has an initial short com- 
ponent followed by a longer component. A few represen- 
tative recoveries are shown in the Fig. ISM-31 We have 
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Table SM-1: Energy of various magnetic state with respect 
to the FM state in meV/Ir. 



Configurations 


FM 


AFMl 


AFM2 


AFM3 


GGA+U+SOC 


0.0 


34.55 


-9.16 


29.20 



Figure SM-4: Three possible magnetic ordering for the HP 
phase. 



fitted the recoveries with the formula 

1 - m{t)/mo = A[Bexp{-t/TL) + (1 - B)exp{-t/Ts)] 

(1) 

where and Ts are the longer and shorter components 
of Ti, respectively, and i? is ^ 0.7. varies from about 
40 s at 100 K to about 400 s at 8K. Below 8K, Tl be- 
comes even longer though we did not take detailed data. 
The major differences of Ti between the AP and HP sam- 
ple are, (i) Ti is almost three orders of magnitude higher 
for the HP sample compared to the AP sample and (ii) 
relaxation rate (1/Ti) increases with decreasing temper- 
ature for the AP sample but decreases with decreasing 
temperature (T < 45 K) for the HP sample. 



II. DETAILS OF THEORETICAL CALCULA- 
TIONS 



Density functional theory calculations have been per- 
formed using three different methods, namely (a) plane 
wave based method as implemented in Vienna ah initio 
simulation package (VASP) P, Q (b) full potential lin- 
earized augmented plane wave method (FP-LAPW) [1 
and (c) N*^ order muffin-tin orbital (NMTO) method 
To begin with, we have checked the relative stabil- 
ity of various magnetic configurations using VASP. The 
calculations were performed within local (spin) density 
approximation (LSDA) with generalized gradient cor- 
rection (GGA) of Perdew-Burke-Ernzerhof [H including 
Hubbard U and spin-orbit coupling SOC. The kinetic en- 
ergy cut off of the plane wave basis was chosen to be 
500 eV. Brillouin-Zone integration have been performed 
using 10 X 10 X 4 k-mesh for the AP phase and 10 x 10 x 10 
k-mesh for the HP phase. In order to calculate the DOS, 
the partial DOS, the orbital and magnetic moment cor- 
responding to the lowest energy state, we have employed 
the all-electron FP-LAPW method. We have taken the 
value of on-site Coulomb repulsion (U) as 3 eV to under- 
stand the effect of correlations. Finally to extract the 
tight binding parameters and to understand the crys- 
tal field splitting, we have constructed low energy model 
Hamiltonian using the downfolding procedure as imple- 
mented in the NMTO method Q. 

We have first investigated whether magnetism can 
be stabilized in the AP and the HP phase using 
GGA+U+SOC calculations in the framework of LAPW 



method. For the AP phase, we have considered four dif- 
ferent magnetic configurations, namely FM (both the in- 
tra and inter dimer couplings are ferromagnetic), AFMl 
(both the intra and inter dimer couplings are antifer- 
romagnetic), AFM2 (intra dimer coupling is ferromag- 
netic and inter dimer coupling is antiferromagnetic), and 
AFM3 (intra dimer coupling is antiferromagnetic and in- 
ter dimer coupling is ferromagnetic). 

Our calculations as summarized in Table ISM- II reveal 
that AFM2 state has the lowest energy. To understand 
the microscopic origin of the insulating behavior (see Fig. 
3(a) of the main paper), we performed NMTO downfold- 
ing calculations and constructed a Ir-d only low-energy 
Hamiltonian by integrating out the high degrees of free- 
dom other than Ir-d for the AP hexagonal phase. Diago- 
nalization of the on-site block give the following eigen 
states with energies -4.26, -4.26, -4.11, -0.48 and 
—0.48 eV respectively. 

|1) = -0.8135|a;?/) + 0.5815|yz) 

|2) = -0.5815|a;z) - 0.8135|a;2 - 2/2) 

|3) = \z^) 

|4) = -0.5815|a;?/) -0.8135|yz) 

|5) = -0.8135|a;z) + 0.5815|x^-?;2) 

Clearly, the first three states form the manifold and 
the other two form the e^ manifold. Because of the distor- 
tion in the octahedra, the three fold degeneracy of the t2g 
block is lifted. The t2g-eg crystal field splitting is quite 
large (~ 3.6 eV) which is also visible in the PDOS (see 
Fig 3(a) of the main paper). As the exchange splitting is 
much smaller compared to the crystal field splitting, Ir 
(d^'^) remain in the low spin state. The origin of this gap 
can be understood by ignoring the spin mixing term of 
the SOC Hamiltonian and only concentrating on the mi- 
nority t2g bands. With this simplification, we can diago- 
nalize the SOC Hamiltonian, having matrix element 
L.S|j), (i,j = |1), |2), |3)) and find that the degeneracy 
of |1) and |2) is lifted upon inclusion of SOC, while |3) 
remains unperturbed. This is schematically illustrated 
for a Ir dimer in the inset of Fig. 3(a) of the main paper. 
Our NMTO downfolding calculation also reveal that the 
hopping strength between various i2g orbitals of the two 
Ir's which form the structural dimer is strong. Therefore 
these orbitals will interact strongly with each other and 
form six non-degenerate states (see inset of Fig. 3(a) 
of the main paper). The two Ir will have total three 
electrons in the minority channel and will occupy three 
states among the six non-degenerate levels leading to an 
insulator as explained in the letter. 

For the HP cubic phase, to keep our computation 
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less expensive, we have simulated Ba2lr(Y3/4lri/4)06 magnetic configurations considered here (see Fig. ISM-4[) 
which is very close to the experimentally synthesized to be degenerate. 
Ba2lr(Y2/3lri/3)06 system. Our calculations find all the 
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